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The aim of the paper is to propose a new type of fits in terms of invariant quantities for finding 
the entries of the CKM matrix from the quark sector, by using the mathematical solution to the 
reconstruction problem of 3 x 3 unitary matrices from experimental data, recently found. The 
necessity of this type of fit comes from the compatibility conditions between the data and the 
theoretical model formalised by the CKM matrix, which imply many strong nonlinear conditions on 
moduli which all have to be satisfied in order to obtain a unitary matrix. 
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I. INTRODUCTION 

In the description of the electroweak interactions 
within the Standard Model a fundamental role is played 
by the unitary CKM matrix, parametrised by four inde- 
pendent parameters: three mixing angles, Oij, and a CP 
violating phase, S. On the other hand from experiments 
one finds numerical values for other quantities such as 
the moduli, or the angles of so called unitarity triangles, 
and by consequence the reconstruction of unitary matri- 
ces from experimental data is a central problem in the 
electroweak interactions. Only recently a mathematical 
solution to this problem was found, fT] , and one aim of the 
paper is to explore its implications upon the electroweak 
phenomenology. One conclusion of the above cited pa- 
per was that it is necessary to provide phenomenological 
models for analysing data in terms of rephasing invari- 
ant quantities, i.e. the CKM matrix moduli, the Jarl- 
skog invariant, J, 0], the area of any unitarity triangle, 
Q, or the phases of its complex entries In this pa- 
per we will choose the moduli as independent parameters 
because with this choice the reconstruction of a unitary 
matrix from data is essentially unique, see Ref. In the 
same time this choice allows us to define phenomenolog- 
ical models which lead to sets of non linear constraints 
which have to be satisfied by data in order to come from, 
or to be compatible with a unitary matrix. Hence the 
data could be compatible with unitarity, or even could 
disprove the CKM unitary model, i.e. nobody guaran- 
tees us that the moduli \Uij\, the mixing angles 9ij, or 
the phase 6 extracted from experiments come from a uni- 
tary matrix. In phenomenological analyses it is usually 
assumed that irrespective how the measured data are, 
they are compatible to the existence of a unitary matrix, 
and in literature one finds statements such as: the most 
stringent test of CKM unitarity is the relation 



ud 



V v. = 1 



(1) 



see, e.g., Ref. @, where Vij are the measured moduli val- 
ues, statement which is not true, because there is a nat- 
ural embedding of unitary matrices into a larger class of 



matrices, that of double stochastic ones, [6|, and we have 
to find methods to discriminate between the two sets. 

A 3 X 3 matrix M is said to be double stochastic if its 
elements satisfy the relations 



> 0, 



E 



(2) 



The unistochastic matrices, are a subset of the dou- 
ble stochastic ones defined by 



1%j — I Uij I 



(3) 
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where t/ is a 3 x 3 unitary matrix. Hence even if all the 
six relations such as ([T]) are exactly satisfied the corre- 
sponding matrix could be non-unitary. 

The double stochastic matrices have an important 
property, they are a convex set, i.e., if M\ and Mi 
are double stochastic, so is their convex combination 
olM\ + (1 — a)M2, a S [0,1], as it is easily checked. 
This property is very important because it was the nec- 
essary ingredient to devise a method for doing statistics 
on unitary matrices. 

The relations ^ together with the embedding rela- 
tions ([3]) define a phenomenological model which allows 
to find formulae for the mixing angles Sij — sin 9ij and 
cos (5 in terms of four independent moduli Vij, and the 
compatibility condition between data and unitarity prop- 
erty leads to the most constraining unitarity condition, 
namely, cos S as function of moduli should take physical 
values, i.e. — 1 < cos 6 < 1, see Ref. [l|. 

The choice of four moduli as independent parameters 
is very appealing from a theoretical point of view, since, 
as we said before, in this case the reconstructed unitary 
matrix is essentially unique. On the other hand it is 
naturally to assume that the physical results of any phe- 
nomenological analysis of experimental data must be in- 
variant with respect to the choice of the unitary matrix 
form, e.g., Kobayashi-Maskawa, 0, or Chau-Keung, [1], 
form. Although both the above cited parameterisations 
depend on mixing angles, 9ij and CP non-conserving 
phase S, only the last parameter is invariant, the numeri- 
cal values for 9ij depend on the chosen form. The mixing 
angles being not invariant quantities, it is better to avoid 
their use in phenomenological analyses. If the physical 
results are not invariant, this could be a signal that the 
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unitary model is not compatible with experimental data, 
or that our phenomenological data analysis is done with 
wrong parameters. 

The above considerations imply that for defining a 
novel type of CKM fits we have to make clear a few nat- 
ural assumptions, which we will state as axioms. Con- 
cerning the experimental side an obvious axiom could be 
the following: 

1) the numerical values for all the measured moduli, \Uij\, 
should be the same irrespective of the physical processes 
used to determine them. 

The second axiom comes from the constraints imposed 
by data on the explicit theoretical tool used for doing 
phenomenological analyses. Thus the second axiom could 
be: 

2) the physical results obtained from data analyses must 
be invariant with respect to the choice of four independent 



invariant quantities used to parametrise the data. 

Because in this paper we will use moduli as invariant 
parameters, and the number of four independent moduli 
groups equals 58, one gets 165 different forms for cos 5, 
which when computed by using real data one finds differ- 
ent values, instead of a single one. Thus the above axiom 
will have in this paper the form: 

2a) the physical results must be invariant with respect to 
all the choices of four independent moduli groups used to 
parametrise the data. 

Of course for defining phenomenological models and 
for a full reconstruction of unitary matrices from error 
affected data we have to use one explicit form of the 
CKM matrix, and in the following we use the form from 
Ref. which, by making full use of its invariance at 
multiplication at left and at right with diagonal phases 
matrices, we write it as 



C12C13 C13S12 Sl3 



U = -C23Si2e** - C12S23S13 ci2C23e** - S12S23S13 S23C13 (4) 



3125236"^ - C12C23S13 -Cl2S23e"^ - S12C23S13 C23C13 



with the standard notation: 



= cos 9ij and Sij = 



sin 9ij; 9ij, ij = 12, 13, 23, denote the mixing angles and 
5 is the phase that encodes the CP-violation in the elec- 
troweak sector. 

The paper content is as follows. In Section 2 we define 
two phenomenological models that will provide the nec- 
essary and sufficient conditions the data have to satisfy 
in order that they should come from a unitary matrix. 
Eventually by finding that data come from a unistochas- 
tic matrix we have to provide an algorithm for the recon- 
struction of U from the error affected data. In Section 
3 we test the moduli data available from experiments, 
lattice computations, and global fits. We show that the 
experimental errors are quite large so there exists a con- 
tinuum of unitary matrices compatible with data. The 
paper ends by Conclusion. 



II. PHENOMENOLOGICAL MODELS 

By phenomenological model we understand a relation- 
ship between the entries of a unitary matrix and the mea- 
sured quantities, such as the moduli. Hence in the fol- 
lowing we assume the knowledge from experiment of the 
moduli of unitary matrix entries such as (jH), which we 
write as the entries of a positive matrix 



V 



I ' ud 'us ' ub 

vl vl 

' cd ' cs ' cb 



(5) 



^td 



^ts ^tb 



For the current state of the art concerning the determina- 
tion of the above quantities see, e.g., Refs.H and [i]. As 
the notation suggests we make a clear distinction between 
the elements of the unitary CKM matrix U and the pos- 
itive entries matrix V provided by data. In other words 
we make a distinction between the theoretical quantities 
Uij and the experimental moduli Vij , although in an ideal 
situation the relation \Uij\ = Vij will hold. 



The main phenomenological problem is to see in what 
conditions from a matrix such as ^ one can reconstruct 
a unitary matrix as (jH), i.e. to see if the data are com- 
patible with the theoretical model. The compatibility 
conditions are better understood in the frame of a phe- 
nomenological model, which in our case will be a rela- 
tionship between the theoretical object ^ and the ex- 
perimental data ([5]). The theoretical tool we have at our 
disposal for defining such a relationship is the unitarity 
property of the CKM matrices. The aim of any phe- 
nomenological analysis is at least twofold: a) checking 
the consistency of data ([5]) with the theoretical model 
Q, and b) determination of theoretical parameters Sy 
and S from the experimental data ([S]) if the data and 
the theoretical model are compatible. In other words our 
scope is the reconstruction of the unitary matrix (|4]) from 
experimental data. In the following we present two phe- 
nomenological models: unitarity condition method and 
unitarity triangle method. 

In both the models we use the embedding of unitary 
matrices into the set of double stochastic matrices, see 
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([3]), implying that the following relations should hold 

i—d,s.b i—d^s,b 

E \Uv\^-^= E K:'-1 = 0, j = d,.s,b (6) 

z— ii,c,t i—u,c,t 

In the above relations T^j are numbers affected by errors 
obtained from experiments, and the exact fulfilment of 
Eqs.® is not at all guaranteed in any phenomenological 
analysis. Even if they are exactly satisfied, the unitarity 
property could not be satisfied. And this is one main rea- 
son for using different notations for experimental moduli 
Vij, and for theoretical moduli \Uij\. The distinction is 
also necessary since the phenomenologists are not aware 
of the natural embedding of unitary matrices, through 
relation into a larger set, that of double stochastic 
matrices ([2]), i.e. they are not yet aware of the necessity 
to find a criterion for the separation of unitary matrices 
within the double stochastic set, only the intersection of 
the two sets being relevant for the electroweak physics 
phenomenology. 

A. Unitarity condition method 

The specific condition defining this model is the fulfil- 
ment of the relation 

V=\Uf 

between the theoretical object (|H) and the experimental 
data ([5]), relation that has to be understood as working 
entry wise leading to the following relations 

1/2 _ „2 2 ^r2 _ 2 2 T^2 _ 2 

^ud ^ '-12'-13' ''us — *12'-13' ''lib — ^13 

T/2 _ 2 2 y2 _ 2 2 

''ch ^ '^23''13' ^tb — '-13'^23' 

= ■S12C23 + S?3S23Cl2 + 2S12S13S23C12C23 COS (5, 

= C12C23 + «12Sl3^23 - 2S12S13S23C12C23 COS (5, (7) 

^td = ■S13C12C23 + S12S23 - 2S12S13S23C12C23 cos (5, 

^ts = S?2Sl3C23+C?2S23 + 2si2Sl3S23Cl2C23COS(5 

This phenomenological model is similar to that pro- 
posed by Wolfenstein, [13], i.e. it is a direct relationship 
between the measured values Vij and the theoretical pa- 
rameters entering the unitary matrix ([4]), and 5. The 
difference between the two approaches is the following: 
we make no approximations on the right hand of Eqs. ([7|. 
and we will explicitly make use of the double stochastic- 
ity relations ([6]). One sees also that on the left hand 
in Eqs. ([7]) there are sets of nine numbers, V^^, obtained 
from experiments, and on the right hand enter only Jour 
independent parameters. Hence the consistency problem 
of the equations ([7]) is a natural one, and it has to be 
resolved. An other remark is that the 5 dependence of 
equations ([7]) is through the cosine function, which is an 
even function, so we can restrict the range of 5 to (0, tt). 



when we use as independent parameters the moduli Vij , 
without loss of generality. 

The main problem to be solved is to find in what con- 
ditions the system ([7]) could have a physical solution, for 
arbitrary numbers Vij satisfying Eqs.®, i.e. from the 
full set of double stochastic matrices. In Ref. [1| it was 
shown that the necessary and sufficient conditions the 
data have to satisfy in order to the matrix ([5]) comes 
from a unitary matrix are 

< < 1, and -\<cos5<\ (8) 

when Sij and cos 5 are found by resolving the equation 
system ([7]). The most constraining condition is the sec- 
ond one, i.e. cosJ G (—1,1), which is the separation 
criterion between the double stochastic and unitary ma- 
trices, whose fulfilment is compulsory. In fact for the 
3x3 matrices, if the relations ([6]) are satisfied, there al- 
ways exists a solution for Sij which is physical, and, e.g., 
it could be obtained by using three independent relations 
from the first five appearing in equations ([7]). The rela- 
tions ([8]) are the consistency conditions between the data 
and the theoretical model. 

To see the constraining power of the above relations 
([5]) we assume for the moment that the relations (O are 
exactly satisfied. The matrix ([4]) depends on Jour inde- 
pendent parameters and we choose them as four inde- 
pendent moduli, i.e. four \Uij\. Now we use the relations 
^ and the entries from the data set ^ to construct a 
matrix, S — (Sij), whose square entries form a double 
stochastic matrix. In order to simplify the formulae form 
we make the notation 

Vud = a, Vus = b, Vub = c 

Vcd = d, Vcs=e, Vcb^J (9) 

For example, if the four independent parameters are the 
moduh Vus = b, Vub = c, Vcd = d, Vcb = /, the matrix 
S"^ = (5-2.), where 

S = iSij)= (10) 

/ VI - 52 - c2 b C \ 

d yi-d^- j^ J 

V Vb^ +C^-d^ Vd2 + f2 _ 52 ^1 _ c2 - /2 / 

is double stochastic. The above choice is motivated by 
the fact that the moduli 6, c and / are the essential pa- 
rameters used by the CKM fitter groups. By using the 
relations V^j = Sfj in equations ^ one gets values for 
Sij and cos S as follows 

5(1) ^ ^ ^ 

43^ - = 7=^' 4'=K. = c (11) 

V 1 - Kb V 1 - c2 

COS (5(1) = 

d2(l ~ c2)2 - 62(1 _ c2) + /2(62 - + {b^ + c^)) 
2&C/V1 - - C2 
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By doing similar computations, but now with the in- stochastic matrix similar to (jlOp . one gets 
dependent parameters a, 6, d, e, entering in a double 

I 



,(2) 



5(1) 

^3 



5(2) 

'23 



~7a 



62 



cos (5. 



(2) 



- &2) _ ^2(1 _ „2) ^ ^2(^2 _ 52 ^ ^2(^2 ^ ^2)) _ {^a^ _ 1) + ^2(1 ^ ^2)) 



2 a 6 VI - a2 - fe2 ^ 1 _ rf2 _ g2 ^^2" 



P + d'^ + e2 - 1 



(12) 



and so on. Now we provide two numerical matrices, built 
as in pI7|) . whose squared moduli are double stochastic. 



Ml 



993907 
106 

2269 
1W 



22691 
lO^i 

24327 
25000 



V17007251 \ 
106 1 



3V453251 
50000 



%/21545351 %/16271665 V998351289149 



(13) 



10= 



106 



The moduli entering the matrices Mi and M2 have 
been chosen such that they should not be too far from 
the recommended moduli values given by the PDG 2006 
data, i.e. within the error corridors, and not too 
far from the central values of the last fit by CKM fitter 
Group [TT|. Their difference is 



0.973907 0.22691 0.004124 
0.2269 0.97308 0.040394 
0.004642 0.040338 0.999175 



Ml - M2 





0.0002 
-0.00595 





-4.7 X 10-5 
0.00114 



(15) 



M2 



973907 
106 

2267 

V112265351 
106 

0.973907 
0.2267 
0.010596 



22691 

9\/7306909 
25000 

V15364455 
10= 



V17007251 
106 



3V453251 
50000 

\/998260569149 
106 



0.22691 
0.973127 
0.039197 



0.004124 \ 
0.04039 
0.999175 / 



i.e. the entry moduli of the matrix ()15|) are much small 
(14) than the current experimental errors on the correspond- 
ing moduli. In fact only the modulus Afi(2, 1) was 
changed by an amount of 2 x 10~^ and this change was 
propagated to get again a double stochastic matrix. 

By using the formulae (fTTj) and (fT2l) one gets from Mi 



= , '''''' . 0.226912, ^ ^^M^ 

V999982992749 lO^ 



'23 ^ 60^/ ^^^^oonno^.n ^ 0-040395, 1 = 1,2,... (16) 



999982992749 

613855254083724801174896537 

cosd'-^ = , w 1.068 

64224997110656250V8008162692156206007509 



The first remark is that Sy and cos 5 take the same val- 
ues when the data come from an exact double stochastic 
matrix, i.e. they are independent of the four indepen- 
dent moduli we choose to parameterise the data. This 
is a consequence of the fact that the double stochasticity 
(unitarity) properties do not change if we interchange the 
columns and/or the rows between themselves, or we use 



the transposed, or the complex conjugated matrix. The 
second remark is that, although the entries of the matrix 
Ml are good from an experimental point of view, they 
are not not good from a theoretical point of view since 
cos S is unphysical. 

From the second matrix, M2, we get 



5 



JO _ 226910 (.) _ v/17007251 (.) _ r 45325r ^_ 

^'"7999982992749' " 10^ ' '^3 - t^uy 999932992749 ' 



1897412766859734762379196101 

coso*^ ' = , « —0.165 

128449994221312500V8008162692156206007509 



From a theoretical point of view the second matrix, 
M2, is compatible with unitarity. Please remark that 
in both the cases the mixing angles take identical 
values although only M2 comes from a unitary matrix. 
Hence the separation criterion between unitary matri- 
ces and purely double stochastic ones is provided by the 
physical values taken by cos 5. By looking at the differ- 
ence matrix Mi — M2 , Eq. p5|) , one sees that the moduli 
differ by numbers which are smaller than the experimen- 
tal errors and, in spite of this, cos 5 has a substantial 
jump from a unphysical value, 1.068, to a physical one, 

- 0.165, which is equivalent to <5 = 99.5°. Thus the com- 
patibility condition between data and theoretical model, 

— 1 < cos (5 < 1, puts very strong conditions on data in 
order to they come from a unitary matrix, and, as the 
above example shows, the condition is very sensitive even 
to small moduli variation. 

When the relations ^ are exactly satisfied the re- 
construction algorithm of a unitary matrix from a dou- 
ble stochastic one is the following: start with a double 
stochastic matrix as M^, or M|, solve the equations ([7]) 
and obtain results as those given by equations (jl6p and 
(fT7)l . If the numerical value for cos (5 satisfies the in- 
equalities ([H]), then with the values for Sij and cos 6 go 
to equation ^ and find the corresponding unitary ma- 
trix. In the above cases only the numerical results from 
M2 are compatible to the existence of a unitary matrix, 
while those from the Mi matrix are not, although the 
matrix Mi is a double stochastic one, and from the usual 
phenomenological point of view both data could be con- 
sidered as being "physical". The reconstruction of the 
unitary matrix from the M2 data is the matrix U which 
is obtained by the substitution of numerical values for 
and cos (5 given by Eqs. P?|) into formula 

The true real case is when the double stochasticity re- 
lations ([6]) are not exactly satisfied, and by consequence 
the double stochastic matrices, built such as the matrix 
S, are different. To see what happens in this case we will 
use only the central values from the PDG 2006 data 0, 
which are 

a = 0.97377, 6 = 0.2257, c = 4.31 x 10"^ (18) 
d = 0.230, e = 0.957, / = 41.6 x 10"^ 

and those from the fit 12], namely 

a = 0.97504, 6 = 0.2221, c = 3.505 x 10~^(19) 
d = 0.222, e = 0.97422, / = 40.8 x lO'^ 



With them we form the double stochastic matrices cor- 
responding to the phenomenological models (jlip and, re- 
spectively, p^ . and one gets in the first case 

cos (5^^) = 25.98, cos 6^"^^ = 1.58 (20) 

and 

cos (5(1) = 0.6, cos (5(2) = -0.377 i (21) 

in the second case. If in the first case the results could be 
considered as being "normal" since the moduli were ob- 
tained by doing some weighted means on moduli data ob- 
tained by experimenters, in the second case the numbers 
were obtained from a fit which used the unitarity trian- 
gle approach, formalism that is supposed to take properly 
into account the unitarity property. Unfortunately this 
does not happens, the second form for cos 6 providing 
an imaginary value. As we will show in the next section 
this is a characteristic of this fit, not an accident. The 
corresponding values for cosi5 are different because, e.g., 
in the first case, the parameters, a, 6, c, d, e and / do 
not come from the same doubly stochastic matrix, e.g., 
/ 7^ Vl — rf^ — e^, or numerically, 0.0416 7^ 0.1768, and 
so on. Thus even in the case when the relations ([6]) are 
exactly satisfied, by choosing four independent moduli 
and constructing with them a double stochastic matrix, 
the numerical results show that unitarity could not be 
satisfied, many cos 6 values being not physical, and these 
values depend on the chosen four independent moduli, al- 
though the solution for cos 5 of Eqs. Q must be unique. 
Hence we have to find a solution to this problem. 

B. Unitarity triangle method 

The second phenomenological model is defined by the 
orthogonality relations of rows, and, respectively columns 
of a unitary matrix, together with the double stochastic- 
ity relations (O. The last condition was never used in 
the previous approaches, see Ref.[l3|, or |l3]. Although 
there are six such relations usually one considers only 
the orthogonality of the first and the third columns of U, 
relation that is written as 

UudK, + U^aU:, + UtdU:, = (22) 

where * denotes the complex conjugation. The above 
relation can be visualised as a triangle in the complex 
plane, and usually it is scaled by dividing through the 
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middle term such that the length of one side is 1. In fact 
one may divide by any other term because all the three 
possible triangles are similar, i.e. they have the same 
angles. On the other hand all the six triangles obtained 
from all orthogonality relations such as (|22p are equiva- 
lent, i.e. they have the same area. A, and furthermore 
the relation J = 2A holds, where J is the Jarlskog in- 
variant, see The side lengths of these triangles are 
auxiliary parameters, and have no physical significance. 
In contradistinction to sides, the angles are measurable 
quantities such that they are important from an exper- 
imental point of view, see in this respect Ref. • By 
taking into account the relation (fTO)) . the other sides of 
the triangle ([22|) have the lengths 



R 



(1) 



UudU*f, 



R 



(1) 



'cb 



dT 



(23) 



df 



and, respectively. 



(2) 



R 



R?' - 



UudU*h 



Ucdu:, 



^^IHH, (24) 

Vl - - rf^ Va^ + 62 -t- c;2 + e2 - 1 
dVl-d^-^ 



Since the exact fulfilment of Eqs.® does not mean the 
fulfilment of the unitarity property, as our numerical re- 
sults ()16lll7p show, that property is implemented in this 
approach by the conditions: all ratios should be posi- 
tive, Ra > 0, a = b, t, . . . , and the following inequalities 
should be satisfied 3 



\R 



R 



< 1< i? 



R. 



(0 



z = l,2. 



(25) 



that are equivalent with the conditions ([5]). Indeed by 
computing the above ratios by using matrices Mi and 
M2 one gets 



R 



(1) _ 



0.438, Rl 



(1) _ 



0.506, 



^ = 1.156, i?f ^ = 0.438 



such that the relations (pS)) give 

0.068 < 1 < 0.944, 0.718 < 1 < 1.594 



(26) 



(27) 



showing the perfect (theoretical) equivalence at this level 
of both the phenomenological models. If now we compute 
the above ratios by using moduli from the fit [l2| one 
finds 



r['^'^ = 0.83, 



Ri^^ = 0.377, 



i?f ^ = 0M3i, Rf^ = 0.404 



(28) 



in perfect agreement with the results (pij) . showing the 
inconsistency of the fit. 



The Eqs. (|23ll24)) represent the correct form for the side 
lengths of the "standard" unitarity triangle. Its angles 
are easily related to the phases of Ucd and Utd- 

The other unitarity triangle which will provide another 
two independent angles is given by the orthogonality of 
the second and the third columns 



One gets similarly 



UusU* 



R' 



(1) 



be 



(29) 



(30) 



f^l-d^~ P 



and, respectively 



UcsU 



cb 



e \/\ ~ (P- — e^^ 



(31) 



VcsUlb 



fe2 + (i2 _|_ g2 _ 2 



Computing the above ratios with the same numbers as 
in ((28l) one gets 



= 0.02, = 1.01, 

i?!^) = 0.038i, = 1.01 



(32) 



showing the same phenomenon as that given by the re- 
lations (|28p concerning the implementation of unitarity 
constraints for this triangle. 

Similar to the first case, two angles of the triangle gen- 
erated by the relation ([29]) are directly connected to the 
phases of the elements Ucs and Uts- From the above re- 
lations we may obtain all the angles of the two triangles 
generated by the relations ((22| . and, respectively, ((29)) . 

For getting the relationship between the phases of 
complex entries of the CKM matrix ([4]) and the mod- 
uli and phase S we make the following notation Uki — 
±|C/fc; je*'^'"' , fc = 2, 3, ^ = 1,2, where the minus sign will 
be taken in front of IC/21I and IC/32I. We make full use 
of the property that in our parameterisation ((4]) of the 
CKM matrix, the entries Uud, Uus, Uub, Ucb and Utb 
are all real quantities, and of the invariance of unitarity 
triangles angles with respect to the scaling factor. Thus 
we write the relations (|23p in a complex form as 



cos W21 + i sin W21 _ Ucd U*f^ 
Rb Uud U*i^ 

(coscc;3i + ismuj3i) Rt _ Utd Utb 
Rb Uud U*. 



(33) 



Please remark that the phases 0^21 and W31 coincide with 
the phase of —Ucd, and, respectively, Utd- The above 
phases coincides, in some cases modulo tt, with the angles 
of the unitarity triangle obtained by scaling the relation 



7 



(|22)) dividing by the first term, i.e. UudU*i^. Since the 
triangles are similar the following relations hold W21 — 
7, W31 — a. From the first relation ([33]) we easily obtain 



sm UJ21 



cos 0^21 



S12 523 C23 sin(5-Rb 

Sl3 C12 

S23(si3 S23 C12 + S12 C23 cos d)Rb 



(34) 



where from one gets 
tanaj2i 



S13 C12 



S12 C23 sin (5 



si3 S23C12 + S12 C23 cos (5 



(35) 



The above formula for tanu;2i depends only on theoreti- 
cal parameters entering and does not depend on the 
lengths of the unitarity triangle. It shows that it is not 
necessary the use of unitarity triangles construction to 
obtain their angles, they can be easily obtained from the 
moduh Vij, through relations such as (fTT|) and (fT2|) . 

From relations (p4|) we get a direct formula for cos 5. 
Indeed, from the identity sin^ 0^21 + cos^ a;2i = 1 we find 



^2 2 
cos5=£i3^ 



S23(s 



23'>''12 '-23 



'13 ■'23 



2S12S13 S23 C12C23 



(36) 



that depends on Sij and Rb- If in it we substitute the 
formulae for Sij from relations pT|) , and Rb from relations 
we find the expression given in Eg. pTjl for cos 5, 
that shows at this level the perfect equivalence between 
the two approaches. Hence there are explicit formulae 
for cos 5 in terms of triangle lengths and mixing angles 

Sij ■ 

Doing similar calculations by starting with the second 
Eg. ([55)1 one gets 



tancj3i : 
and, respectively, 



S12 S23 sini5 



-si3 C12 C23 + S12 S23 cos 6 



'23 



'13 '-23 '^12) 



2 S12 Si3 S23 C12 C^3 



(37) 



(38) 



In conclusion the orthogonality relation ([221) provides 
only two independent angles, that depend on four in- 
dependent parameters, S12, S13, S23, and 5, and for a 
complete determination of parameters entering a unitary 
matrix we have to use another unitarity triangle, for ex- 
ample that generated by the relation ([^ . By doing sim- 
ilar computations one gets 



tanu;22 



tan LJ32 



C12 C23 sin 5 



-S12 si3 S23 + C12 C23 cos (5 ' 
C12 S23 sin (5 



S12 S13C23 + C12 S23 cosi5 
Hence we can write the matrix (U)) under de form 



U 



Uud 



Uus 

-\UtsW' 



u. 



ub 



(39) 
(40) 

(41) 



and we call u>ij the fundamental phases since with them 
one gets the angles of all unitarity triangles. The phases 
entering the matrix (j4T|) corresponding to the moduli ma- 
trix M2 are given by 

W21 = 99.46°, W31 = 121.427°, (42) 
W22 = 99.50°, W32 = 98.153° (43) 

Our proposal is to use the matrix form (|4ip in all phe- 
nomenological analyses, and the experimenters have to 
measure all its phases which really appear as parameters 
in many electroweak processes. For a similar proposal 
seeRef.i. 

The important conclusions of this subsection are: a) for 
a complete determination of a unitary matrix from data 
when one uses the unitarity triangle model one needs 
the use of at least two unitarity triangles if one wants 
to obtain reliable results, b) the use of double stochas- 
ticity relations ^ is compulsory, otherwise the fit gets 
senseless, and c) checking of fit results to see if them are 
invariant with respect to the choice of four independent 
moduli used to parametrise the data. The numerical re- 
sults, provided by the relations (pij) . (^5]) and ([5^ . show 
that all the above requirements were not fulfilled by the 
fit fl3. 



Recovery of unitary matrices from error 
affected data 



Until now we supposed that the double stochasticity 
relations ([6]) were exactly satisfied, i.e. the data were not 
affected by errors. Even in this case by using measured, or 
fit determined moduli, and by forming with them double 
stochastic matrices we got that the numerical computa- 
tions on cos (5 and lengths of unitarity triangles lead to 
unphysical results, see Eqs. pl^ . ^(MM) and By 
consequence we have to see how the reconstruction al- 
gorithm which works for data coming from exact double 
stochastic matrices has to be modified in order to provide 
reliable results in the presence of errors. The experimen- 
tal data on moduli do not satisfy the double stochastic- 
ity relations ^ , hence the first condition which must be 
imposed when doing a fit is that these relations should 
be satisfied with a great precision, otherwise one cannot 
speak of unitarity fulfilment, because the recovery pro- 
cess gets senseless. Since the double stochastic matrices 
obtained from data by using different groups of four in- 
dependent moduli are really different, leading to different 
values for cos (5, see e.g. Eas. (j20ll2T]) . we have to impose 
that all cos 5 values should be (approximately) the same. 
On the other hand the explicit form for cos 5 depends on 
the four independent moduli we choose to parameterise 
the data, hence in analysing the data we have to make 
full use of the axiom 2a) from Introduction. 

The recovery method of unitary matrices from data 
that we expose in the following is a least squares method 
for checking the compatibility of data with the theoreti- 
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cal models in both the approaches, and stresses the ne- 
cessity that the x^-function have to contain two kind of 
terms: the first has to impose the fulfilment of unitarity 
constraints, and the second should take into account the 
physical quantities measured in experiments. The piece 
of unitarity constraints has an independent part provided 
by the double stochasticity relations, and a dependent 
one upon the phenomenological model. Constraints im- 
plied by double stochasticity relations are written as 



Xc/s 



(44) 



E E y^-A + E E 



j—u,c^t \i—d,s,b 



j—d^s,h 



The constraints generated by the unitarity condition 
method are given by 

= ^(C0S<5W -C0Sp))2, -1<C0S(5«<1 (45) 
i<j 

A similar formula one gets for the second phenomenolog- 
ical model where instead of cos 5 one uses the side lengths 
Ra, see Ref.[l|. The third component which takes into 
account the experimental data has the form 



xl 



E 



(46) 



is a double stochastic matrix. In other words we have 
to do statistics on the set of double stochastic matrices 
generated by unitary ones. Thus this property allows 
us to calculate correctly mean values, (M), and error 
matrices, (Tj\/, for a set of double stochastic ones, and by 
consequence for unitary matrices, as follows 



(M) 




(49) 



/ k—7i 



J2 \Uk\^]/n {My 



\k=l 



where the square roots are taken entry wise. If the mean 
value matrix obtained in this way is not too far from 
a unitary one, one can reconstruct from it an (approxi- 
mate) unitary matrix. Also important is the relationship 
between the "central" values matrix and those obtained 
by adding the error matrix ctm- Thus the above formulae 
suggest that the ±3 a matrices should be calculated by 
the formulae 



M+ = \JaP + 3cr2^, ^ _ 2,alj (50) 

In this way we provided a complete formalism for doing 
fits on data from the electroweak sector, which is a very 
robust one such as can be seen in the next Section. 



where di are the theoretical functions one wants to be 
found from fit, dt is the numerical matrix that describes 
the corresponding experimental data, while a is the ma- 
trix of errors associated to di. Hence a x^-test could be 
the function 



X 



xL 



2 I 2 
Xl +X2 



(47) 



which will be used in numerical computations. 

The compatibihty condition of Eqs. ([7]), together with 
the invariance of physical quantities with respect of the 
four independent moduli group chosen to parametrise the 
data, imply strong conditions on all the moduli given by 
the relations coscJ'*-* ~ cosJ^^^ for i ^ j, see Eg. (H5t . 
which are not easy to implement in a fit, but gives us at 
least one reward: a method for doing statistics on unitary 
matrices. 

It is well known that the problem of doing statistics on 
(moduli) of unitary matrices was an open problem, see 
e.g. [ll]. The embedding relations ^ suggest that the 
right quantities to look at them are the square moduli, 
I C/ij P . Indeed the convexity property of double stochastic 
matrices implies that if we have a set of unitary matrices 
Ui, . . . ,Un then 

i—n i=n 

= ^ x,- |C/fcp,, = 1, (48) 

1=1 1=1 
< a^i < 1, i — i, ■ ■ ■ ,n 



III. TESTING DATA AND FITS 

In the following we want to show what are the main 
consequences of the above formalism, especially what are 
the subtleties and the results obtained when doing statis- 
tics on unitary matrices by the above method. For that 
we will make use of all 165 different forms for cos 5. We 
test the Reviews of Particle Physics 2004 ^ and 2006 
data 0, the lattice computations of CKM moduli p7j . 
and the published fits [III and [l^ . 

1. 2004 PDG data [3]. As it is weU known there 
one finds upper and lower bounds for each modulus, at 
95% confidence level. We used the mean values, calcu- 
lated as half the sum of upper and lower bounds, and the 
corresponding cr,; computed as half the difference of the 
bounds. One gets that the double stochasticity is quite 
well satisfied by the mean values, all the six sums differ- 
ing from 1 by amounts of order 10~^ — 10"'*, which lead 
to 



(cos (5) = 0.063 -t- O.OOSi 



(51) 



Because the mean real and imaginary parts are small, the 
above results suggest that the central data are consistent 
with a 5 value around 90°, and a fit done by using our 
method provided the numerical matrix 



0.9748379 
0.222699 
0.0098077 



0.222885 
0.974025 
0.0399759 



0.00365667 \ 
0.0409987 
0.999153 j 



(52) 



9 



which leads to the following value 

6 = (90.01 ±0.41)° 



(53) 



For the above matrix the double stochasticity property 
is well satisfied, all six relations ([6]) taking values whose 
magnitude is of the order of 10~^, and this implies a small 
(statistical) error for S. 

If one uses the recommended values [lB| for the entries 
on the first two rows, and the mean values for moduli on 
the third row one gets big non-physical values 



(cos J) = -24.15 - 8.16i, cTc 



504.21 - 0.39i (54) 



which show that the recommended values are quite far 
from values compatible with unitarity requirements. 

2. 2006 PDG data We made a data modification 
for one single modulus, Vtb, by taking its central value 



equal to Vtb = 0.99912, instead of Vtb = 0.771^:24, since 
we saw that the obtained results are similar, and one gets 



(cos (5) = 21.557- 21.320i, CTc 



62.738 + 7.326i (55) 



which shows that the mean data are far from being com- 
patible with unitarity. A fit with our formula (|47p around 
the central values from PDG 2006 data provides the ma- 
trix 



0.974290 0.225260 0.004114 \ 
V-,, = I 0.225064 0.973579 0.038602 
0.010256 0.037441 0.999246 / 



where from we get 



(101. 18± 0.18)" 



(56) 



(57) 



The 2006 PDG data ^9] bring an innovation: entering 
of global fits results into Review of Particle Physics book, 
see Ceccucci et al. contribution to CKM quark-mixing 
matrix 1 1811 in the book. The fits taken into account are 
[m and |l9l|. By using the central values provided there 
one finds 



(cos (5) = 0.5084, CTc 



0.1002 



which is equivalent to 



S = (59.44: 



-6.46\c 
-6.93) 



(58) 



(59) 



which shows at the level of la an important statisti- 
cal spreading caused by a poor fulfilment of the double 
stochasticity relations for the first row and the first col- 
umn. If we look at the central values Vc + cr, and respec- 
tively Vc — (T, where +a and —a are slightly different, see 
Eq.(11.26) in 1^, we get 



(cos 6^ 



1.4175 + 0.0508i, 
6.1236 - 0.0118i 



and respectively 

(cos (5^) 

^cos 6^ 



0.3076 
6.2063 



0.2324i, 
0.0115i 



(60) 
(61) 

(62) 
(63) 



Both results show that moduli values 14 ± ct have no 
physical relevance and provide one numerical example 
showing the shortcomings of the unitarity triangle ap- 
proach such as it is used in the present day. The big 



values for (t„ 



5 5± 



confirm that many cos S values take 



values outside the physical range (—1, 1), and an impor- 
tant number takes imaginary values. For example in the 
second case, Vc — cr, the real values are within the interval 
cos (5 e (—16.45, 19.17), and the imaginary ones, 19 val- 
ues from 165, within the interval cos (5 G (—0.355, 7.525)i 

In our opinion this example shows clearly the two 
causes which generate such results: a) fitting only with a 
single unitarity triangle and using essentially only three 
moduli, 5, c and /, and b) non existence of a sound pro- 
cedure for doing statistics on moduli of unitary matrices. 

3. By using the lattice results from Ref. JJ^ on CKM 
matrix moduli one gets 

(cos(5) = 8.962 - 11.315i, acosS = 20.15 + 5.03i (64) 

which shows that lattice computations are still far from 
results compatible with unitarity constraints. The above 
results show that the double stochasticity property is not 
satisfied, hence the numerical effort has to be done to 
improve it. By looking for a compatible unitary matrix 
around the mean moduli values from Ref. IitI, with our 
method one gets the numerical matrix 



0.974330 0.225086 0.0041314 \ 
V3= I 0.225003 0.973562 0.0393922 
0.007401 0.038911 0.999215 / 



(65) 



which leads to S = (56.04 ± 0.62)° . This moduli matrix 
can be used as a good "unitary witness" to improve the 
numerical algorithms used to obtain results as those given 
in Ref. flf\. 

The sensitivity of our method to small moduli varia- 
tions can be seen by looking at the entries of matrices 
([56|) and ([65]) . The difference of the moduli on the first 
two rows is of the order of 10~^, while S changes from 
101° to 56°. 

4. CKM fitter Group 2001 results [13]. The real and 
imaginary values for cos S are within the bounds cos S S 
(-1.082, 1.596), respectively, cos (5 e i(-2.431, 0.835) and 
one gets 



(cos (5) = 0.373- 0.087 i, 



0.358 + 0.907i (66) 



Hence one could say that the central values together with 
their ±cr companions are not compatible with unitarity. 
This happened because the invariance of the physical re- 
sults with respect to all the choices for the four indepen- 
dent moduli groups, and the double stochasticity rela- 
tions were not yet implemented in the unitarity triangle 
model, and, by consequence, the fit have no physical sig- 
nificance. 

5. CKM fitter Group 2005 results [Il|. The last pub- 
lished results by CKM fitter Group are, in our opinion, 
the best published results until now. Thus our fourth 
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selected matrix is given by the central values from [llj 
see their Table 3 



0.97400 0.2265 0.00387 \ 
Va = I 0.2264 0.97317 0.04113 

0.00826 0.04047 0.999146 / 



(67) 



The errors on moduli are of the order 5.5 x 10~^ for Vud 
and Vcs, and of order 2.5 x 10"'^ for Vus and Vcd- The 
double stochasticity is well satisfied obtaining numbers 
of the order (10~^ — 10~®). On this matrix one can check 
how the fulfilment degree of double stochasticity prop- 
erty reflects in cos 5 values. One gets for 5 and Jarlskog 
invariant, J, the following values 



5 = {&2mtlZ) 



J = (3.84t°:^^) X 10 



(68) 



We remind that for a double stochastic matrix the sta- 
tistical error of 5 is zero, such that the errors for 5 give 
a measure of the double stochasticity fulfilment. As a 
curiosity we remark that, from the point of view of uni- 
tarity constraints, the numerical results (|67j) are better 
than those appearing in PDG 2006 data, [l^, perhaps 
since the last ones were obtained by merging results from 
two different fitting groups. Unfortunately the errors ac- 
companying the matrix (|67p have no physical relevance 
if are used according with the present day rules. 

By doing a fit with our method around the values from 
Eq. ([S7)) one gets such a matrix, which is only slightly 
different from the above one. If we check now the results 
for V4 -f Icr, and V4 — lu, where the matrices +cr and —a 
are different, see their Table 3, and one adds the errors 
such as they are given, one gets 



(cos (5)+ = -2.214 - 4.559i, 
(cos(5)_ = -2.91- 5. 33i 



(69) 



If we proceed by using our method, i.e. we calcu- 
late 5± by using the moduli obtained with the formulae 
\/V^ ± one gets 



(61.02+^;^")° and 5^ = (61.91 



+3.45NC 
3.56^ 



(70) 



Also one finds that the -|-cr error matrix is overestimated 
si nce, for ex ample, computations with the moduli matrix 
\/Vil + 3(7^ lead to complex values for (cos (5), namely 
(cos(S) = 0.554 - 0.015i. 

6. Other selected matrices. Around the central val- 
ues from PDG 2006 data we imposed a supplementary 
constraint, namely, Vtd/Vts = 0.208, and we got 



0.973556 
0.228289 
0.00851678 



0.228407 
0.972704 
0.0409463 



0.00431 
0.0416 
0.999125 



(71) 



and from it one finds 5 = (63.706 ± 0.228)°. 

More exotic matrices are: a) around central PDG 2004 
data values 



0.975282 
0.220702 
0.0106931 



0.220925 
0.974397 
0.0418558 



0.00340635 \ 
0.04299 
0.999066 / 



(72) 



with 

5= (98.1 ±6.4)° 
b) around PDG 2006 data central values 



0.97328 0.2295 0.0062 \ 
T/7 = I 0.229112 0.97256 0.04044 (73) 
0.0153 0.037939 0.99916 / 



which leads to 



5 ^ (155.7 ± 13.5)" 



or c) 



equivalent to 



0.973514 0.228585 0.004311 \ 
0.228564 0.97264 0.0416 (74) 
0.005316 0.041484 0.999125 / 



5 = (6.4 ±4.0)' 



The double stochasticity property is not very well satis- 
fied by the last three exotic matrices. We provided them 
to see how our method for doing statistics on unitary ma- 
trices works. However, from a phenomenological point of 
view, all the above eight matrices are good, and each 
one leads to a perfectly acceptable unitary matrix. The 
above examples show that the obtained unitary matrices 
depend strongly on the numerical values around one looks 
for a compatible unitary matrix. Changing them a little 
bit one gets different values for 6. Hence the data are 
not known with a sufficient precision in order to tighten 
significantly the CP violating phase 5. 

The novelty brought by the embedding of unitary ma- 
trices into the convex set of double stochastic matrices is 
the following: we can use the double stochasticity prop- 
erty, and from the (approximate) double stochastic ma- 
trices Vi, i = 1 . . . , 9, where Vg denotes the approximate 
numerical form of the M2 matrix, we get a continuum of 
(approximate) double stochastic matrices 



(75) 



< a;, < 1, i = l, 



1=1 



< 1 



continuum which depends on eight arbitrary parame- 
ters Xi. Its consequence is that with the W matrix one 
can obtain practically any value for S within the inter- 
val (0°, 180°), all of them being relatively on the same 
footing in what concerns their agreement with the exper- 
imental data on the moduli Vij , and with the fulfilment 
of unitarity constraints. 
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If we use the formulae (|49|) to do statistics on the 
above nine matrices we find at the symmetric point 
Xi — 1/9, i — \, - ■ ■ ,8, the matrix 

/ 0.974134 0.225921, 0.0042982 \ 
W = 0.225766 0.973323 0.0409179 (76) 
\ 0.009583 0.040019 0.999153 / 

and its associated error matrix 




0.00119 0.00118 7.4 x 10"^ \ 
0.00118 0.00117 1.0 X 10--* (77) 
5.7 X 10-5 11.68 X 10-5 9.88 x 10"^ j 



which provide the foUowing value for 5 

(5 = (81.12 ± 1.23)° (78) 

Although the moduli matrices, 1^, lead to 5 values prac- 
tically on all the physical interval, (0, 180)°, the mean 
value W satisfies the double stochasticity property quite 
well shown by the relatively small statistical error in for- 
mula (TfT]) . which leads to a small error for 5, see Eg. ([78| . 
Of course the true matrix a has to be calculated by taking 
into account all Ui matrices, by using the same property: 
one takes the mean of squared af, which lead to a new 
matrix whose entries are a little bit bigger than those 
appearing in ([77]) . 

The above results show that the present day fits are 
not able to determine the CKM matrix entries with a 
high precision, and by consequence neither 5, nor the J 
invariant can be well defined. The computations, as well 
as all the numbers obtained by the so called global fits, 
use only partial information coming from experiments, 
which has as consequence the existence of a continuum 
of unitary matrices compatible with a set of experimental 
data and unitarity constraints. As the numerical compu- 
tations show this type of global fit gives better results and 
in the same time it provides a sound method for doing 
statistics on unitary matrices. 

IV. CONCLUSION 

In the paper we have shown that the compatibility be- 
tween the moduli data and the theoretical model ([4]) can 



be obtained if and only if three conditions are simultane- 
ously satisfied: a) the double stochasticity relations ([6|) 
hold to a level of at least 10"^, b) the physical condition 
— 1 < cos (5 < 1 is satisfied, and c) the numerical values 
for cos 5 do not depend on the four independent moduli 
chosen to parametrise the data. 



Our numerical checking of the two approaches in Sec- 
tion 2 has shown that both the theoretical models, unitar- 
ity triangle method and the unitarity condition method, 
send the same signal concerning the compatibility of data 
with unitarity requirements. Thus now we have two phe- 
nomenological approaches, and both of them have to give 
results consistent each other. 

The numerical implementation of the x^-test, such as 
(j47|) . is a little bit more complicated in both the ap- 
proaches when one makes use of the exact form of uni- 
tarity constraints. This happens since by looking at the 
denominators of the relations (fTlIT2ll . ((23lfT7)) . p0ll3T|) . 
and those similar to them, one sees that there appear 
square roots which by modifying a little bit the mod- 
uli values, change easily their values from real to imagi- 
nary ones. On the other hand by using approximations, 
as those used in the standard form of unitarity triangle 
model, one easily finds non physical results. However, 
as our numerical computations from Section 3 show, the 
implementation of all the exact unitarity constraints can 
be done. 



Our computations have shown that the phenomenol- 
ogists are not yet aware of the constraining power of 
unitarity, this one requiring a moduli matching to an 
order of 10"^ — 10~®, much lower than the experimen- 
tal errors, tuning that can not be obtained by using ap- 
proximate formulae. The influence of errors on the final 
results is that in case they are not small enough, the dou- 
ble stochasticity property provides us continuum sets of 
approximate unitary matrices. Combining these results 
with a correct method for doing statistics on unitary ma- 
trices we have obtained a powerful tool for checking uni- 
tarity properties of data that will become available at the 
LHC machine in the next years. 
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